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ABSTRACT 

Mid-term quasi-periodicities in solar coronal mass ejections (CMEs) during the most 
recent solar maximum cycle 23 are reported here for the first time using the four- 
year data (February 5, 1999 to February 10, 2003) of the Large Angle Spectrometric 
Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO). In 
parallel, mid-term quasi-periodicities in solar X-ray flares (class >M5.0) from the 
Geosynchronous Operational Environment Satellites (GOES) and in daily averages of 
Ap index for geomagnetic disturbances from the World Data Center (WDC) at the 
International Association for Geomagnetism and Aeronomy (IA GA ) are also examined 
for the same four-year time span. By Fourier power spectral analyses, the CME data 
appears to contain significant power peaks at periods of ~ 358 ± 38, ~ 272 ± 26, 
~ 196 ± 13 days and so forth, while except for the ~ 259 ± 24-day period, X-ray solar 
flares of class > M5.0 show the familiar Rieger-type quasi-periods at ~ 157 ± 11, 
~ 122 ±5, ~ 98 ± 3 days and shorter ones until ~ 34 ± 0.5 days. In the data of daily 
averages of Ap index, the two significant peaks at periods ~ 273 ± 26 and ~ 187 ± 12 
days (the latter is most prominent) could imply that CMEs (periods at ~ 272 ±26 and 
~ 196 ± 13 days) may be proportionally correlated with quasi-periodic geomagnetic 
storm disturbances; at the speculative level, the ~ 138 ± 6-day period might imply 
that X-ray flares of class > M5.0 (period at ~ 157 ± 11 days) may drive certain types 
of geomagnetic disturbances; and the ~ 28 ± 0.2-day periodicity is most likely caused 
by recurrent high-speed solar winds at the Earth's magnetosphere. For the same three 
data sets, we further perform Morlet wavelet analysis to derive period-time contours 
and identify wavelet power peaks and timescales at the 99 percent confidence level for 
comparisons. Several conceptual aspects of possible equatorially trapped Rossby-type 
waves at and beneath the solar photosphere are discussed. 

Key words: oscillations — space weather — Sun: activities — corona — coronal 
mass ejections — magnetic fields 



1 INTRODUCTION 

Mid-term quasi-periodicities (one to several months or 
longer) in various diagnostics of solar flare activities and 
sunspot numbers or areas etc. during a few years around the 
solar maximum phase have been extensively searched for and 
monitored at many electromagnetic wavelengths (Rieger et 
al. 1984; Kiplinger et al. 1984; Dennis 1985; Ichimoto et al. 
1985; Delache et al. 1985; Bogart & Bai 1985; Bai & Sturrock 
1987; Oliver et al. 1988; Ribes et al. 1987; Lean & Brueckner 
1989; Ozgiie & Atac 1989; Lean 1990; Carbonell & Ballester 
1990; Droge et al. 1990; Pap et al. 1990; Kile & Oliver 1991; 
Verma et al. 1992; Ballester et al. 1999; Cane, Richardson, & 
von Rosenvinge 1998). These activities of observational re- 



search were triggered by the landmark discovery by Rieger 
et al. (1984) of a ~ 154-day quasi-period in solar 7-ray flare 
rates registered by the Gamma-Ray Spectrometer (GRS) on- 
board the Solar Maximum Mission (SMM) two decades ago. 
Besides this oft-quoted quasi-period of ~150-160 days, there 
are other notable quasi-periods around ~128, ~ 102, ~ 78, 
and ~ 51 days during the maxima of different solar cycles 
(Dennis 1985; Bai 1992; Bai & Sturrock 1991) from various 
data sets. Empirically, these quasi-periods seem to hint at 
quasi-subharmonics of equatorial solar rotation period (e.g., 
Sturrock & Bai 1992; Bai & Sturrock 1993). 

As to the plausible physical origin of such quasi- 
subharmonics with periods longer than the solar rotation 
period, we proposed recently (Lou 2000a, b) that, in the 
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statistical sense, the existence of equatorially trapped large- 
scale (scales comparable to Rq) Rossby-type waves (Rossby 
ct al. 1939; Papaloizou & Pringle 1978; Provost et al. 1981; 
Saio 1982; Wolff & Blizard 1986; Wolff 1998) may be ulti- 
mately responsible for quasi-periodically modulating or trig- 
gering smaller-scale magnetic catastrophes in solar active 
regions. In this scenario, solar magnetic active regions be- 
come, recurrently, vulnerable to catastrophes via statisti- 
cal accumulations of magnetic stresses and energies through 
magneto-convective turbulence. 

In terms of specific observational diagnostics, we sug- 
gested (Lou 2000b) that the photospheric magnetic flux 
emergence might be triggered via magnetohydrodynamic 
(MHD) buoyance instabilities (Parker 1955) modulated by 
equatorially trapped Rossby-type waves owing to the pres- 
ence of large-scale subphotospheric magnetic flux (Gilman 
1969; Lou 1987; see Ballester et al. 2002 for most recent 
emerging magnetic flux observations), that large-scale co- 
herent velocity patterns or "cells" may exist over the so- 
lar surface (see Beck, Duvall, & Scherrer 1998 and Ulrich 
2001), and that more precise measurements of solar sur- 
face elevation may reveal < 0.1" variations (see Kuhn et 
al. 2000 for much smaller 100-m high "hills" with an az- 
imuthal separation of [8.7 ±0.6] x 10 4 km ). As Rossby-type 
waves involve fairly slow and large-scale vortical oscillatory 
disturbances with periods longer than the solar rotation pe- 
riod, we speculated (Lou 2000a, b) that an uninterrupted 
time sequence of large-scale coronal mass ejection (CME) 
events might contain similar mid-term quasi-periodicities. 
The LASCO/SOHO global observations of solar CMEs in 
space offer an unprecedented and unique opportunity for 
such an investigation. Uninterrupted LASCO/SOHO CME 
data is now available for four years (February 5, 1999 to 
February 10, 2003). 

Meanwhile, we examine possible mid-term quasi- 
periodicities for powerful X-ray solar flares of class >M5.0 
and daily averages of Ap index for geomagnetic storm dis- 
turbances for several reasons. First, as a proxy of 7— ray 
flares (Rieger et al. 1984), strong X-ray flares are known to 
contain Rieger-type periodicities during past solar maxima. 
We would like to reconfirm this feature and provide a test of 
our current data analysis at the same time. Secondly, strong 
X-ray flares might correlate with global CMEs in some sta- 
tistical sense. Thirdly, in terms of both space weather and 
solar terrestrial interactions, we would like to establish pos- 
sible correlations among CMEs, intense X-ray flares, and 
geomagnetic disturbances characterized by daily averages of 
Ap index. Systematic investigations on quasi-periodicities 
may reveal such correlations in a unique manner. 

For all three data sets, we shall first use Fourier anal- 
ysis to identify possible quasi-periodicities. To complement 
the Fourier spectra, we further carry out Morlet wavelet 
analyses to derive contour plots in period-time domains for 
the same three time series of CME, X-ray solar flares, and 
geomagnetic Ap index data, respectively. With a few ex- 
ceptions, major peaks identified in the wavelet analyses are 
roughly consistent with major peaks in the Fourier power 
spectral analyses. Additional information of when certain 
quasi-periodicities occur or recur in different data sets may 
provide valuable clues for possible physical connections. 

In broader astrophysical contexts, we note in passing 
that global MHD tidal waves including Alfven-Rossby type 



Figure 1. (a) Daily counts of CME events from LASCO for the 
four-year time span shown on top. (b) Fourier power spectrum 
of the CME data, (c) Comparison of statistical probability for 
power larger than a preset value between the CME data sequence 
(asterisk *) and an artificially yet sufficiently randomized CME 
data sequence (+). 



waves (Lou 1987, 2000b, 2001) may give rise to valuable pe- 
riodic or quasi-periodic signatures from magnetized compact 
stars such as neutron stars and white dwarfs that are cov- 
ered by corotating thin yet dense "plasma oceans" involving 
extremely strong magnetic fields ranging from 10 7 G to 10 15 
G. It is our desire to use the Sun as an astrophysical lab- 
oratory to search for, to identify, and to understand stellar 
Rossby-type waves. 



2 DATA ANALYSIS AND RESULTS 
ON QUASI-PERIODICITIES 

The time sequence of daily CME counts was derived 
from a preliminary list of LASCO/SOHO compiled by 
the Naval Research Laboratory (NRL) at the web- 
site lasco-www.nrl.navy.mil/cmelist.html. The solar X-ray 
flare data were taken from the GOES at the website 
www.ngdc.noaa.gov/stp/SOLAR/solar.html and from the 
website www.sec.noaa.gov/Data/solar.html. The data of the 
daily averaged Ap index for geomagnetic storm disturbances 
were obtained from the WDC of the I AG A at the website 
swdcdb.kugi. kyot o-u . ac.jp/wdc/Sec3. html. 

2.1 Fourier Power Spectral Analyses 

The uninterrupted data for daily counts of CMEs from 
LASCO/SOHO are displayed in Fig. la. Prior to Feb. 5, 
1999, there were several unfortunate interruptions of SOHO 
operations and those earlier CME data are therefore ex- 
cluded. Fig. lb is the Fourier power spectrum of the data 
shown in Fig. la. In spite of the dominance (Fig. la and lb), 
peak A at ~1101-day period is not reliable as this period is 
not much less than the total data length of 1467 days (edge 
effects). While being preliminary, the next three significant 
peaks B, C, and D at quasi-periods of ~ 358±38, ~ 272±26, 
and ~ 196 ± 13 days, respectively, are longer than the usual 
Rieger-type quasi-periods (Rieger et al. 1984; Dennis 1985; 
Sturrock & Bai 1992) in the flare-related data (see Table 
1 for significant CME power peaks identified). To estimate 
possible errors in periods of identified power peaks, we take 
half of the full width at half maximum (FWHM) about a 
power peak in Fig. lb as the angular frequency variation 
Alu. Then the period variation is given by AP = — PAw/w. 
That is, P ± AP would give an estimated period variation 
range for an identified power peak. In the same manner, we 
have obtained estimates for power spectra of X-ray flares 
(class >M5.0) and of daily averages of Ap index (see Table 
1 for details). By Table 1, AP increases with increasing P 
in general. 

For time series of the daily CME data shown in Fig. 1, 
we took the following steps to estimate the statistical signif- 
icance of Fourier spectral power peaks (e.g., Delache et al. 
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Table 1. List of quasi-pcriodicities (Data sets from February 5, 1999 to February 10, 2003) 



Item 




A 


Identified Quasi- 
B C 


-Periods by Spectral Power Peaks in 
D E F G 


Unit of Days 
H I 


J 


K 


L 


CME 


P 

AP 


1101.22 
671.48 


358.33 
37.68 


271.99 
26.22 


195.88 
13.38 


110.80 
3.26 


100.00 
2.88 


66.25 
1.31 


60.64 
0.90 


57.26 
1.00 


35.85 
0.40 


33.49 
0.39 


20.62 
0.17 


Flares 


P 

AP 


1505.00 
589.46 


259.48 
24.23 


156.77 
10.89 


122.19 
4.88 


98.15 
3.25 


67.69 
2.18 


63.73 
2.34 


42.16 
0.78 


38.71 
0.55 


33.53 
0.52 






Ap 


P 

AP 


940.63 
249.85 


364.11 
38.91 


272.81 
25.55 


187.34 
12.05 


138.07 
6.33 


91.49 
5.19 


60.64 
2.38 


27.80 
0.22 


27.08 
0.21 


24.69 
0.17 


24.00 
0.25 


23.29 
0.14 



All original data were taken from February 5, 1999 to February 10, 2003. Significant power peaks based on Fourier analyses are 
identified and summarized here. We have periodicity= P ± AP and AP = —PAlu/ui where Au> is a half of the full width at half 
maximum (FWHM) about a power peak. See Figs. 1, 2, and 3 in the text for CME, X-ray solar flares, and Ap index, respectively. 



Figure 2. The same panel arrangement as in Fig. 1 but for the 
X-ray solar flare data (class >M5.0) from GOES during the same 
time span indicated on top. 



Figure 3. The same panel arrangement as in Fig. 1 but for the 
daily averaged Ap index of geomagnetic disturbances from the 
WDC/IAGA during the same time span indicated on top. 



1985). First, the daily CME data is picked at random within 
the time sequence and rearranged several to more than ten 
times to form an artificially randomized sequence (similar 
to the process of shuffling cards) with much reduced coher- 
ent periodic signals if any. For such a white noise sequence 
of independent Gaussian distributions with variance a 2 , a 
Fourier transform would yield a power spectrum character- 
ized by a probability distribution p(w) — exp[— w/(2a 2 )] for 
spectral powers greater than a preset value w. The approx- 
imate straight line of symbol + in Fig. lc represents the 
natural logarithm of p(w), \n[p(w)], versus w for an artifi- 
cially randomized CME data sequence with a variance of 
a 2 ^ 1.46 x 1CT 3 {a 3.8 x 1CT 2 ) estimated from the slope 
of this line. The relevant a-levels of fluctuation amplitudes 
are then shown (viz., the two horizontal dashed lines in Fig. 
lb corresponding to 3-tr and 4-a levels as well as a sequence 
of numeral marks on top of Fig. lc). Next, the natural loga- 
rithm of probability for spectral powers greater than a preset 
value w is estimated from the power spectrum of the actual 
daily CME counts (symbol * in Fig. lc). A comparison of 
the two curves in Fig. lc clearly indicates that for powers 
greater than the 3cr-level, the CME data starts to depart 
from random noises. Power peaks B through K in Fig. lc 
are regarded as significant. Periods of peaks J and K are 
~ 36 ± 0.4 and ~ 33 ± 0.4 days, respectively. We ran the 
same test procedure for the other two data sets and ob- 
tained qualitatively similar properties. That is, for both Ap 
index and X-ray solar flares of class > M5.0 data sets exam- 
ined below, deviations from random noises are also notable 
for power peaks > 3cr-level (see Figs. 2c and 3c). 

It should be noted that such a procedure can be exe- 
cuted whatever the count distribution of the original data 
sequence may be. In fact, we have obtained data count dis- 



tributions about the respective means for all three data se- 
quences (i.e., CMEs, X-ray flares, and Ap index). It turns 
out that both CME and X-ray flare data more or less follow a 
Poisson distribution (that is, P^{v) — exp(— jx)^ jv\ where 
/x is the mean count rate and v is the actual count rate; see 
e.g., Taylor 1981), while the Ap data deviate significantly 
away from a Poisson distribution (not surprisingly) . 

Entirely parallel to the panel arrangement of Fig. 1, 
Figure 2 presents a similar analysis on daily counts of X-ray 
solar flares of class > M5.0 from the GOES. By Fig. 2c, the 
estimated variance is a 2 = 3.52 x 10~ 5 (a = 5.9 x 10~ 3 ). As a 
confirmation, the X-ray flare data in Fig. 2a does contain the 
familiar Rieger or Rieger-type quasi-periodicities < 157±11 
days (see Fig. 2b and Table 1 for power peaks identified 
in the solar X-ray flare data) which do not have significant 
counterparts in the CME power spectrum except for periods 
of ~ 34 ± 0.5, ~ 39 ± 0.6, and ~ 98 ± 3 days (Fig. lb and 
Table 1). As we deal with the most powerful X-ray flares, 
this appears to be consistent, in trend, with the fact that 
Rieger-type quasi-periodicities were first discovered (Rieger 
et al. 1984) in 7-ray flares by the SMM/GRS followed by 
detections in soft X-ray flares (Kiplinger et al. 1984). Note 
that peak B of ~ 259 ± 24-day period appears significant 
enough (> 4a level) and may be related to the ~ 272 ± 
26-day period in the CME data. For reasons noted earlier, 
peak A at 1505 ± 589-day period again might not be reliable 
even though its apparent confidence level is merely slightly 
< 95%. 

Similarly, Figure 3 contains the relevant information of 
the daily averaged Ap index for geomagnetic storm distur- 
bances, with an estimated variance a 2 = 0.08 (a = 0.28). 
Power peaks C and D in Fig. 3b at periods ~ 273 ± 26 and 
~ 187 ± 12 days may correspond to the two peaks at peri- 
ods ~ 272 ± 26 and ~ 196 ± 13 days of the CME data in 
Fig. lb. There also exists a fairly prominent peak H (> 4- 
a level) at a ~ 28 ± 0.2-day period. The ~ 138 ± 6-day 
period might correspond to the ~ 157 ± 11-day period in 
X-ray flare data; admittedly, this correspondence is some- 
what weak given the ranges of error estimates. Should this 
correlation be real, then there might exist subclasses of ge- 
omagnetic disturbances which are proportionally correlated 
with major X-ray solar flares (class >M5.0). As already sus- 
pected, the 941 ± 250-day period may not be reliable due 
to edge effects. In sharp contrast, by a visual inspection of 
Fig 3a and 3b, the most dominant power peak D at period 
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Figure 4. Period-time contours of Morlct wavelet analysis on 
the solar CME data (Feb 5, 1999 to Feb 10, 2003) with uiq = 6 
and the normalized power. Timescales (periods) of power peaks 
A through E are summarized in Table 2. The cone of influence 
(COI) is beneath the dashed curve. 



Figure 5. Period-time contours of Morlet wavelet analysis on 
the X-ray solar flare data (class > M5.0) with u>o = 6 and the 
normalized power. Timescales (periods) of power peaks A through 
F are summarized in Table 2. The cone of influence (COI) is 
beneath the dashed curve. 

of ~ 187 ± 12 days for Ap index should be physically real. 
In Fig 3, the prominent peak H around ~ 28 ± 0.2-day pe- 
riod in the daily averaged Ap index is physically identified 
with the solar rotation that recurrently brings high-speed 
solar wind streams from low-latitude coronal holes towards 
the Earth's magnetosphere. Note that this ~ 28 ± 0.2-day 
period is almost absent in the CME data. 

2.2 Considerations of Statistical Significance 

Figs, lc, 2c, and 3c provide a sense of statistical significance 
for power peaks identified in the periodograms of CME, X- 
ray flare, and daily averaged Ap index data, respectively. 
Basically, when a spectral amplitude (square root of a power 
peak) becomes higher than the 3<r-level, the statisitics of 
corresponding power peaks in the periodograms (asterisks *) 
starts to deviate from power peaks derived from randomized 
data series (crosses +). The higher the spectral amplitude, 
the larger the deviation in statistics. 

Any one of the time series analyzed here contains ran- 
dom fluctuations that may not necessarily obey Gaussian 
distribution. The corresponding power spectrum naturally 
contains random variations. For independent daily counts of 
sufficiently large number in a data sequence, the correspond- 
ing Fourier amplitudes obey Gaussian distribution by the 
central limit theorem (e.g., Fan & Bardeen 1995). Therefore, 
the power probability distribution of Z = P(ui) being in the 
interval z < Z < z + dz is Pr{z < Z < z + dz} = exp(— z)dz 
where P(co) is normalized by a 2 . The cumulative distribu- 
tion function is Pr{Z < z} — f* cxp(—z')dz' — 1 — exp(— z). 
Consequently, the probability for Z > z is Pr{Z > z} = 
exp(— z). To search for the maximum value of power peaks 
among N independent frequencies such that Z are indepen- 
dent random variables, the corresponding probability is 

Pr{Z max > z} = 1 - [1 - exp(-z)f . 

By this expression, it follows that a preset threshold z = zo 
for power peaks among N independent variables is related to 
the probability (false alarm probability po) that a peak, be- 
ing produced randomly by chance, is greater than zo (Lomb 
1976; Scargle 1982; Ballester et al. 2002). 

To further assess the statistical significance of power 
peaks in Figs. 1 to 3 in this approach, we estimate the de- 
Figure 6. Period-time contours of Morlet wavelet analysis on the 
Ap index data with u)q = 6 and the normalized power. Timescales 
(periods) of power peaks A through D are summarized in Table 
2. The cone of influence (COI) is beneath the dashed curve. 



Table 2. List of power peaks in Morlet wavelet contours 



Item 


Periods in 


Days for Contour Peaks 






A 


B 


C 


D 


E 


F 


CMEs 


343.0 


187.0 


102.0 


38.3 


36.1 




X-ray Flares 


242.5 


144.2 


72.1 


66.1 


39.3 


25.5 


Ap index 


288.4 


187.0 


66.1 


51.0 







In normalized Morlct wavelet period-time contours of Fig. 4 
(CME data), Fig. 5 (X-ray flare data), and Fig. 6 (Ap index 
data), we identify power peaks and timescales (periods) at the 99 
percent confidence level (1 percent significance), respectively. 



tection threshold zq = — ln[l — (1 — Po) lllN ] in spectral power 
with a false alarm probability po, where N is the number of 
independent frequencies over which a power peak is searched 
for (Scargle 1982). In our case, N = 734 (a half of the total 
number of days of observations) ; for a false alarm probability 
po — 0.01 (99% confidence), we have a detection threshold 
of z = 11.2cr 2 and for p = 0.05 (95% confidence), we have 
a detection threshold of zo = 9.6a 2 . In Fig. 1, peaks B — D, 
F, G, J, K are above the 99% confidence level; and peaks 
E, H, and / are above the 95% confidence level. In Fig. 2, 
peaks B — E and H — J are above the 99% confidence level. 
In Fig. 3, peaks C, D, H, I, and J are above the 99%; and 
peaks B, F, G, and K are above the 95% confidence level. 

The tentative correlation found between mid-term 
quasi-periodicities of the CME and Ap data (Figs. 1 and 3) 
can be of physical significance. Besides the cause of recur- 
rent high-speed solar winds with southward interplanetary 
magnetic fields for geomagnetic disturbances, a geomagnetic 
storm may set in when a CME or a part of a CME impinges 
upon the Earth's magnetosphere (Cane et al. 2000; Wang 
et al. 2002). Among all solar CMEs, there is certainly a 
significant fraction that miss the Earth's magnetosphere ow- 
ing to their various initial onset orientations over the solar 
surface. Statistically, the probability, roughly proportional 
to the solid angle subtended by the Earth's magnetosphere 
towards the Sun, that CMEs hit the Earth's magnetosphere 
is higher for more frequent occurrence of CMEs. It is in this 
proportional sense that periodicities in CME and Ap data 
should correlate with each other. This is an important per- 
spective for probing solar-terrestrial interactions and space 
weather conditions. This empirical correlation also strength- 
ens our confidence that the two low-frequency power peaks 
C and D in the CME data are likely to be physically real. 

The solar X-ray flare ~ 259 ± 24-day period is close to 
the CME ~ 272 ± 26-day period and the Ap ~ 273 ± 26-day 
period. If not coincidental by chance, this might indicate 
a certain causal relation among energetic X-ray solar flares 
( > M5.0) and CMEs. Meanwhile, we caution that it would 
be premature to claim the existence of different types of 
solar flares associated with global CMEs. 



1 For example, Jupiter's hectometric radio emissions and extreme 
ultraviolet aurorae were found to be triggered by arrivals of inter- 
planetary shocks according to simultaneous observations of the 
Cassini and Galileo spacecraft (Gurnett et al. 2002). 
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2.3 Morlet Wavelet Analyses 

As already noted earlier, mid-term quasi-periodicities of var- 
ious solar flare related diagnostics may change for maxima of 
different solar cycles. It is thus of considerable interest to see 
whether such periodicities change within a few years around 
the maximum phase of a solar cycle. It is also important to 
check whether power peaks of Fourier spectral periods and of 
wavelet contour timescales have reasonable correspondence. 
For correlation studies of different time series data, it is cru- 
cial to know at what times certain periodicities occur or 
recur. 

For these purposes, we perform wavelet analysis on the 
three time series for the CMEs, solar X-ray flares, and Ap 
index using the Morlet wavelet function 

Mv) =<K- 1 ' i e^ on e- ri2/ ' 2 (1) 

where ujq = 6 is chosen. The Fourier time period r and the 
Morlet wavelet timescale s are related by 

T = 47rs/Vo + (2 + ^o) 1/2 ] (2) 

(see Table 1 of Torrence & Compo 1998). For our choice 
of luq = 6, it follows from equation (2) that r = 1.033s, 
namely, Fourier periods and Morlet wavelet timescales are 
nearly equal to each other. The wavelet transform suffers 
from edge effects at both ends of the time series. This gives 
rise to a cone of influence (COI) as indicated by regions 
below the dashed lines in Figures 4, 5, and 6. Caused by 
padding zeroes at the beginning and at the end of data se- 
ries, these edge effects usually lead to a power reduction 
within the COI. For further detailed technical information 
of Morlet wavelet analysis, one may visit the website at 
:/ /paos. colorado.edu/research/wavelets/wavelet2. html 
At the 99 percent confidence level (1 percent signifi- 
cance) and using the Morlet wavelet function defined by 
equation (1), contour plots in period-time (or equivalently, 
scale-time) domains with color coding are shown for the 
CME, X-ray solar flares, and Ap index data in Figures 4, 
5, and 6, respectively, and corresponding contour peaks are 
summarized in Table 2. 

Comparing Table 1 and Table 2, we find that CME 
Fourier periods B, D, F, and J at ~ 358 ± 38, ~ 196 ± 13, 
~ 100 ± 3, and ~ 36 ± 0.4 days roughly correspond to CME 
Morlet periods A, B, C, and E at ~ 343, ~ 187, ~ 102, and 
~ 36 days, respectively; X-ray flare Fourier periods B, C, F, 
and / at ~ 259±24, ~ 157±11, ~ 68±2, and ~ 39±0.6 days 
roughly correspond to X-ray flare Morlet periods A, B, C, 
and E at ~ 243, ~ 144, ~ 66, and ~ 39 days, respectively; 
and Ap index Fourier periods C, D, and G at ~ 273 ± 26, 
~ 187 ± 12, and ~ 61 ± 2 days roughly correspond to Ap 
index Morlet periods A, B, and C at ~ 288, ~ 187, and 
~ 66 days, respectively. 

The most prominent match is the ~ 187-day period in 
both CME and Ap index data (this quasi-periodicity is also 
present in corresponding Fourier power spectra). As already 
noted earlier, we suspected that this Ap index periodicity 
was driven by periodic CMEs. These comparisons suggest 
that CMEs and Ap index for geomagnetic disturbances may 
proportionally correlate with each other in both periods and 
time. We categorically emphasize that power peaks in Ap 
index data around 181 — 187 days are so prominent in Figs. 
3a, 3b, and 6, and that they should be physically real for 



whatever origin. Here, the evidence lends support for the 
CME-driven scenario in the statistical sense. 

Wavelet power peak in C— D range of Fig. 6 (Ap index) 
and wavelet power peak in C— D of Fig. 5 (X-ray flares) have 
comparable timescales of ~ 66 days. They occur around the 
same time and may be interpreted as flare-driven periodic- 
ities in Ap index data. There is also a match around peri- 
odicity of ~ 39 days in both X-ray flare and Ap index data. 
However, there is a time difference for the two relevant power 
peaks in CME and X-ray flare data. From such information 
alone, correlations between solar X-ray flares and CMEs at 
these timescales are not immediately clear. 



2.4 A Test of Incomplete X-Ray Flare Data 

Such CME mid-term quasi-periodicities (present in both 
Fourier spectrum and Morlet wavelet analyses) , if confirmed 
by further independent observations, offer novel diagnos- 
tics for probing and understanding the physical origin of 
CMEs (e.g., Hundhausen 1999; Low 2001). Up to this point, 
one major question arising from our comparative data anal- 
ysis is that there are mid-term quasi-periodicities in the 
LASCO/SOHO data for CMEs but these periodicities, ex- 
cept for a few, are largely independent of the Rieger-type 
quasi-periodicities detected in various phenomena associ- 
ated with energetic X-ray solar flares. 

While the SOHO/LASCO data of CMEs including all 
halo events is the most complete one, the X-ray solar flare 
data of the GOES suffer incompleteness because one cannot 
see the other half of the Sun at any given moment. Might this 
be the cause of differences in some quasi-periodicities seen in 
CMEs and solar X-ray flares? We designed a simple test to 
address this issue with the assumption that a strong X-ray 
solar flare occurs at one footpoint of a CME. In this scenario, 
a CME may be related to one flare at most. For example, we 
may observe a CME without seeing an intense X-ray flare 
perhaps because one footpoint of the CME is located on 
the frontside but another is located on the backside of the 
Sun and the corresponding flare happens at the footpoint on 
the backside. By all combinations, M CMEs per day may 
correspond to M + 1 (0,1, ...,M) possible X-ray flare counts 
on the frontside of the Sun. Based on the CME data, we 
then rely on a random number generator to select one of 
the possibilities as observable flare events per day. The key 
features of the resulting power spectrum thus obtained do 
not change significantly for arbitrary trials. This test seems 
to indicate that should a more complete data set of solar 
flares be available, the major quasi-periodic features in the 
power spectrum would remain. 



3 MID-TERM QUASI-PERIODICITIES 

On the basis of our data analyses reported here, it appears 
that, with high probability, mid-term quasi-periodicities do 
exist in association with diagnostics of various solar activ- 
ities (CMEs and X-ray solar flares) as well as geomagnetic 
disturbances (Ap index). Some of the overlapping quasi- 
periodicities might reveal underlying physical causes. We 
shall discuss qualitatively plausible physical connections. 
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3.1 Quasi-Periodicities in Flux Emergence 

Regarding the theoretical suggestion of equatorially trapped 
Rossby waves (Lou 2000a, b) for Rieger-type periodicities 
in flare-related diagnostics, we note that intuitively, solar 
flare rate and emergence of complex magnetic regions of 
sunspots or sunspot groups should relate to each other in a 
statistically correlated manner (Oliver et al. 1998; Ballester 
ct al. 1999, 2002). Empirically, nearly all solar flares are 
found close to sunspots or sunspot groups and often seen as 
brightenings of the pre-existing plages; the observed quasi- 
periodicities (Bai 1992; Oliver et al. 1998) in sunspot areas 
or number of sunspot groups do indeed correlate with those 
of flare occurrence rates. Is this merely a coincidence? Or, 
does this actually hint at an underlying global mechanism? 
We here discuss two seemingly independent aspects, namely, 
the emergence of magnetic fluxes and the occurrence of solar 
flares, even though the two aspects should in some sense be 
related to each other empirically. 

Whether magnetic fields are generated deep in the Sun's 
radiative interior (e.g., Gough & Mclntyre 1998) or at the 
bottom of the solar convection zone (e.g., Rosner & Weiss 
1985), some of these magnetic fluxes must somehow float 
upward through the convection zone and eventually break 
through the thin photosphere by magnetic buoyancy (Parker 
1955, 1979) to buckle upwards above the solar surface, form- 
ing sunspot pairs or groups. In the absence of a global 
quasi-periodic modulation mechanism, such process of mag- 
netic flux emergence should be completely random. The ob- 
served quasi-periodicities in magnetic flux emergence would 
be consistent with (though by no means necessary) the pres- 
ence of large-scale Rossby-type waves of comparable peri- 
odic timescales. Generally speaking, equatorially trapped 
Rossby-type waves (Lou 2000a, b) can dynamically affect 
subsurface as well as emerged large-scale magnetic fields 
(e.g., magnetic active regions) through squeezing, twisting, 
stretching and vortical motions, even though spatial scales 
of (emerged) individual sunspots or sunspot groups are rel- 
atively small. 2 More specifically, such equatorially trapped 
Rossby-type waves may tip off vulnerable regions for the 
onset of magnetic buoyancy instabilities via dynamic cou- 
plings. This plausibly links quasi-periodicities detected in 
sunspot areas or number of sunspot groups (Lean 1990; 
Oliver et al. 1998; Ballester et al. 1999, 2002) with equa- 
torially trapped Rossby-type waves. 



3.2 Quasi-Periodicities in Flare Diagnostics 

The occurrence of solar flares involves sudden releases or 
bursts of considerable magnetic energies accumulated over 
a certain period of time (e.g., Parker 1979). It is almost im- 
possible to predict the location and time of a solar flare. 3 As 
magnetic fields associated with sunspots or sunspot groups 
are strong in strengths and complicated in structures, it is 



2 The process leading to the formation of such intense magnetic 
structures is nonlinear and complicated (Parker 1955). 

3 This situation is analogous to hurricanes, earthquakes, and tor- 
nadoes in that vulnerable environments for their occurrences can 
be identified with relative ease while their specific occurrence lo- 
cation and time are very difficult to pinpoint or predict. 



quite natural to somehow have a significant amount of mag- 
netic energy stored in special ways as magnetic fields gradu- 
ally evolve in time (Low & Lou 1990; Lou 1992). In contrast 
to nanoflares (Parker 1994) which frequently release mag- 
netic energies on smaller scales (e.g., X-ray bright points) 
for heating the lower corona, this accumulation (or storage) 
stage of magnetic energies is necessary for flare phenom- 
ena in general otherwise solar flares would not be so vio- 
lent and explosive (e.g., Zirin 1988). At critical moments, 
small disturbances of whatever origin may trigger mag- 
netic avalanches, leading to productions of energetic par- 
ticles and wide-spectrum electromagnetic radiations. Only 
in this sense, a sustained passage of Rossby-type waves can 
increase statistically the chance of flare occurrence in pre- 
existing vulnerable magnetic complexes and thus impose 
quasi-periodicities to various diagnostics associated with X- 
ray solar flares. 

3.3 Dynamic Feedback Cycle 

Having argued heuristically about possible roles of equatori- 
ally trapped Rossby-type waves in triggering the emergence 
of magnetic fluxes and the onset of solar flares in a sta- 
tistically correlated manner, we need to address the press- 
ing question of Rossby-type wave excitation (Lou 2000a, b). 
Whether such waves (with much weaker magnitudes) ex- 
ist or not during the solar minimum phase is not presently 
known. Nevertheless, preceding discussions require their 
presence as a large-scale coordinating mechanism during the 
solar maximum phase. 

As a solar flare or a CME goes off, a fraction of their en- 
ergies will back react on the solar photosphere and appear in 
the form of traveling disturbances. Collectively, a fraction of 
total flare or CME energies may provide a necessary energy 
source for generating and sustaining equatorially trapped 
Rossby-type waves during a few years around solar maxima. 
In this scenario, a dynamic feedback cycle may be established 
to sustain quasi-periodicities in magnetic flux emergence and 
flare activities. Failure of sustaining such a feedback cycle 
leads only to occasional periodicities in flare activities (see 
Fig. 5). This may explain why during maxima of some solar 
cycles, the 150-day quasi-periodicity becomes so prominent 
(Oliver et al. 1998) while in others it could be almost ab- 
sent with yet different periodicities more overwhelming (cf. 
Bai 1992; Bai & Sturrock 1991, 1993; Sturrock & Bai 1992) 
- the dominant periodicity can be different during solar 
maxima of different solar cycles (e.g., Bai & Sturrock 1991; 
Bai 1992) and for different diagnostics (Oliver et al. 1998; 
Ballester et al. 1999, 2002). Fig. 5 clearly shows time vari- 
ations of periodicities even within a few years of the recent 
solar maximum phase (cycle 23). This proposed feedback 
scenario leaves considerable rooms for possible correlations 
among intensities/rates of solar flares, selection of periodic- 
ities, and distribution of sunspot groups on both northern 
and southern sides of the solar equator during solar maxima 
and therefore also requires further explorations. 

3.4 Quasi-Periodicities in CME Diagnostics 

CMEs from the solar limb have been routinely monitored 
by various coronagraph observations (cf. Hundhausen 1999 
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for SMM results). With the ongoing Solar and Heliospheric 
Observatory (SOHO) mission in space, even relatively weak 
CME halo events 4 can now be readily detected as well. By 
combining SOHO and ground-based observations, it is pos- 
sible to construct a fairly complete time sequence for CMEs 
occurrence. 

CMEs represent a major class of large-scale solar activ- 
ities involving magnetic field and perhaps sub-photospheric 
magnetic flux emergence. Physically, CMEs are character- 
ized by large-scale eruptive mass losses from the Sun into the 
solar wind (Hundhausen 1999) and are thought to be caused 
by a sudden loss of magnetic equilibria (Low 1990) as pho- 
tospheric magnetic field structures evolve gradually in time. 
The primary energetics of CMEs is believed to be magnetic 
in nature. At a critical stage of large-scale magnetic struc- 
ture evolution, a CME may be triggered by disturbances. 
As a systematic and quasi-periodic source of large-scale dis- 
turbances at the photospheric level, equatorially trapped 
Rossby-type waves may modulate CMEs in a periodic yet 
statistical manner. 



4 APPLICATION OF ROSSBY-TYPE WAVES 

Our comparative observational research is motivated by the 
known mid-term quasi-periodicities, namely, Rieger-type pe- 
riodicities, in various solar flare related diagnostics and our 
analysis here on X-ray flares (class > M5.0) during the re- 
cent solar maximum of cycle 23 confirms once again the exis- 
tence of such Rieger-type periodicities. What is then the un- 
derlying physics for such "discrete quasi-periods" arranged 
roughly like subharmonics and significantly longer than the 
solar rotation period if they are indeed real? Solely from 
the perspective of period matchings, it appears that solar 
Rossby-type waves are the best bet for such discrete subhar- 
monic quasi-periodicities, although many questions remain 
to be addressed for this proposal (Lou 2000a, b). 

Now given the global nature as well as randomness of 
CME events, how does one infer the underlying physics for 
mid-term quasi-periodicities (see Fig. 1 and Table 1) longer 
than those of the Rieger type? The reality of such mid-term 
quasi-periodicities for CMEs appears to be indirectly sup- 
ported by the presence of similar quasi-periodicities in the 
daily average of Ap index (see Fig. 3 and Table 1). The 
basic rationale for such a remote connection is that quasi- 
periodic CMEs should lead to quasi-periodic geomagnetic 
disturbances when a fraction of CMEs encounter the Earth's 
magnetosphere. It is certainly possible to match the iden- 
tified quasi-periods of CMEs with periods of Rossby-type 
waves (Lou 2000a, b) as we did for Rieger-type periodic- 
ities. However, as shown in Table 1, longer periods have 
larger error ranges comparable to the equatorial solar rota- 
tion period, that is, the accuracy is insufficient to warrant a 
one-to-one correspondence without ambiguities. 

With these qualifications and limitations in mind, we 
shall nevertheless discuss several empirical, intuitive and 
conceptual aspects of mid-term quasi-periodicities in refer- 
ence to equatorially trapped Rossby-type waves. 

4 The strengths of these halo CMEs are not necessarily weak. 
They are not as easy to detect merely because they originate 
somewhere within the solar disk. 



For Rieger-type periodicities, we may estimate the 
quasi-periods of equatorially trapped Rossby-type waves by 
taking m = 12, 10, 8, 6 and 4 with n = 1 or 2 in expressions 
(13) and (15) of Lou (2000b). We suspect that the spa- 
tial distribution of sunspot groups or active regions along 
the equatorial zone might be one important factor in the 
selection of the two integers m (k x = m/U® is the az- 
imuthal wavenumber) and n (number of nodes for a Rossby 
wave function along a longitude) during the solar maximum 
phase. By definition, a solar minimum is characterized by a 
fewer number of sunspots and thus less activities. During a 
solar maximum in contrast, radio images of the Very Large 
Array (VLA) as well as X-ray images from Yohkoh reveal 
two quasi parallel chains of 5 or 6 emission-intense "blobs" 
from sunspot groups or magnetic complexes on both north- 
ern and southern hemispheres. 5 As these "blobs" are sites 
of frequent flare occurrence, excited equatorially trapped 
Rossby-type waves are likely to concentrate their powers in 
the relevant m which roughly corresponds to the number of 
"blobs" around the Sun. For ~5 or 6 active centers in the 
solar disk facing us, m would then be ~10 or 12 around the 
circumference. In turn, Rossby-type wave disturbances will 
statistically trigger more flares in active regions. In essence, 
this is part of the dynamic feedback cycle outlined earlier. 

Likewise, the gross symmetry for the numbers of active 
regions across the equator would correspond to n — 2. Thus, 
equatorially trapped Rossby-type waves with n = 2 might be 
favorably excited given the presence of two belts of sunspot 
groups across the equator. It should be noted that for wave 
modes with even n, the north and south are symmetric with 
respect to the equator; and for wave modes with odd n, 
the north and south are antisymmetric with respect to the 
equator. Of course, for a mixture of Rossby-type wave modes 
with both odd and even n, it is not possible to identify a 
symmetry with respect to the equator in the strict sense. 

Observationally, the distribution of active regions is 
sometimes not grossly symmetric with respect to the equator 
during a certain epoch of a solar maximum phase and this 
situation may correspond to odd values of n. This might be 
the cause of north/south asymmetries sometimes observed in 
periodicities of flare rates (e.g., Bai 1987). One should keep 
in mind that for n = 1 or 2 and m > 4, the angular wave pat- 
tern speeds uj p = uj/m of equatorially trapped Rossby-type 
waves are approximately proportional to ~ 2f2 /m 2 . The 
larger the value of m, the slower the angular wave pattern 
speed lj p . By including effects of surface elevation (see equa- 
tions 1-3 of Lou 2000), periods of Rossby-type wave crest 
bumping into the next adjacent active region are 27r/|mw p | 
as specifically given by expressions (13) and (15) of Lou 
(2000b), namely, 

P r Si P Q {\m\/2 + (2n + l)fi 0J R Q /[|m|( 5 .D) 1/2 ]} (3) 
for Rossby-wave periods, and 

P r _ p £ {\m\ + [m 2 + 8n Q Re/(gD) 1/2 ] 1/2 }P Q /4 (4) 
for Rossby- Poincare- wave periods, respectively, where Pq = 

5 In some cases, the separation of different "blobs" may involve 
certain ambiguities. Sometimes, one might hesitate as whether to 
count a "blob" at the limb or not. For the present purpose, we 
mainly focus on the basic concept of such mode selection. 
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2-k/Q.q is the Sun's sidereal rotation period, QqRq ~ 
2 km s _1 is the equatorial solar rotation speed and (gD) 1 ^ 2 
is the solar surface gravity wave speed with g ~ 2.7 x 
10 4 cm s -2 being the solar surface gravity and D (of a few 
hundred kilometers) being an effective thickness of the pho- 
tospheric layer. 

The solar photosphere is magnetized with well-known 
inhomogeneities on smaller scales such as concentrations 
of intense magnetic fibrils (~ 10 3 G) along boundaries of 
supergranules whose typical diameters are on the order of 
~ 3 x 10 9 cm and appearances of sunspots or sunspot groups 
(of sizes comparable to that of a typical supergranule) with 
magnetic field strengths up to 3 — 4 x 10 3 G. Subsurface hor- 
izontal magnetic fields Bh (e.g., Babcock & Babcock 1955) 
may well be patchy on large scales of the order of L ~ Rq 
and might be fairly strong deep in the solar interior (e.g., 
Gough & Mclntyre 1998). However, as the gas density p 
increases rapidly with depth by gravitational stratification, 
the value of the Alfven wave speed Ca = -Bh/(47rp) 1/ ^ 2 may 
not be very large (say, Ca < a few km s _1 in the overshoot 
layer suspected to be strongly magnetized). 

In terms of large-scale magnetohydrodynamic (MHD) 
wave propagations in a thin layer with a relatively small Ca 
and a relatively large azimuthal perturbation scale, the role 
of subsurface magnetic field Bh may be crudely assessed 
(Lou 1987). The presence of Bh introduces a large-scale 
horizontal magnetic pressure force in addition to the quasi- 
hydrostatic pressure force pgr] caused by the surface eleva- 
tion n and thus effectively increases the surface wave speed in 
the form of cl ~ (gD + C\) 1 ^ 2 where D < 500km is an esti- 
mated photospheric layer thickness and g = 2.7 x 10 4 cm s~ 2 
is the solar surface gravity. Replacing the surface wave 
speed c = (gD) 1 / 2 by cl here, expressions (7) and (11) of 
Lou (2000b) then give increased frequencies of equatorially 
trapped Kelvin- and Poincare-waves (thus shorter periods 
than a couple of days), whereas expressions (13) and (15) 
of Lou (2000b) (or equations (3) and (4) here) show that 
the frequencies of equatorially trapped Rossby-type waves 
remain more or less unchanged because they are primarily 
determined by the solar rotation rate Qq and the minor cor- 
rection term involves the speed ratio Q,QRQ/(gD + Ca) 1 ^ 2 
which is reduced further by the presence of Bh . For Rossby- 
type waves, this allows more nearby frequencies to be packed 
into the frequency groups labeled by m with different values 
of n. 



5 OTHER CIRCUMSTANTIAL EVIDENCE 
FOR ROSSBY-TYPE WAVES 

While solar Rossby-type waves or r-modes have been pro- 
posed along several lines of research (Papaloizou & Pringle 
1978; Wolff 1998; Lou 1987, 2000a, b), there is no direct 
solid evidence so far for their detection at the solar pho- 
tosphere except for a few tantalizing cases discussed below. 
We intend to stimulate further investigations of this problem 
both observationally and theoretically . 

The velocity correlation analysis on the data from the 
Solar Oscillation Investigation /Michelson Doppler Imager 
(SOI/MDI) revealed stationary "long-lived velocity cells" 
(Beck et al. 1998) located along the solar equatorial zone. In 
terms of horizontal surface velocities (a few meters per sec- 



ond) and spatial scales (~ 50° in longitude), it is inevitable 
on theoretical ground that such large-scale "velocity cells" 
should travel in the form of Rossby-type waves relative to 
the Sun at speeds slower than the solar rotation (Lou 1987, 
2000a, b). An azimuthal scale of ~ 50° may correspond to 
an azimuthal wave integer m = 6 or 7 which would give an 
equatorially trapped Rossby wave period > 77 but < 102 
days (Lou 2000b). The diagnostic approach of Beck et al. 
(1998) is promising, although more sophisticated velocity 
correlation analyses of the SOI/MDI data are required in 
order to extract signals of large-scale, slowly drifting Rossby- 
types wave patterns over the solar photosphere. 

Using the solar limb data from the SOI/MDI, Kuhn et 
al. (2000) reported periodic signals of 100 m high "hills" that 
are separated by an azimuthal scale of (8.7 ± 0.6) x 10 4 km. 
Note that this surface elevation of ~ 100 m is some 700 
times smaller than the upper limit of 0.1" estimated by Lou 
(2000b). Uncertainties in such an estimate involves the ef- 
fective thickness D of the solar photosphere and the extent 
of partial cancellations in the horizontal velocity divergence 
Vi • v±. The more precise result of Kuhn et al. (2000) casts 
serious doubts on earlier attempted measurements of varia- 
tions in the solar diameter (e.g., Delache et al. 1985; Ribes et 
al. 1987; Yoshizawa 1999). If such solar limb elevations are 
indeed induced by Rossby waves with an azimuthal wave in- 
teger m = 50, then the corresponding equatorially trapped 
Rossby wave drifts at an extremely slow speed relative to 
the Sun with a period of ~ 625 days. It is somewhat sur- 
prising that limb elevations caused by Rossby-type waves 
with smaller values of m were not detected by the same 
technique of Kuhn et al (2000). 

Combining the data from the SOI/MDI and the Global 
Oscillation Network Group (GONG), Ulrich (2001) obtained 
persistent patterns of torsional oscillations with azimuthal 
structures characterized by 1 < m < 8 yet without detect- 
ing signals of quasi-periodicities. For this range of m values, 
periods of equatorially trapped Rossby waves should lie in 
the range of ~ 27 — 102 days (Lou 2000b). It is a challenge 
through a helioseismological data analysis to identify un- 
ambiguously signatures of solar Rossby waves that involve 
large-scale and long-lived structures of flow velocity and sur- 
face elevation over the slowly rotating solar photosphere. 

Ballester et al. (2002) recently examined several his- 
torical records of solar photospheric magnetic fluxes. These 
records include (1) the Mount Wilson total magnetic flux 
(MWTF) from 1966 to 2000 for total daily magnetic flux 
measured in units of 10 22 Mx, (2) the Kitt Peak magnetic 
flux (KPMF) from 1975 to 2000 in units of 10 22 Mx in re- 
gions with magnetic field strength greater than 25 G within 
the latitude band S70°-N70°, and (3) the daily magnetic 
plage strength index (MPSI) and the Mount Wilson sunspot 
index (MWSI) from 1970 to 2000. They derived both pe- 
riodogram and wavelet contours. The results of their data 
analysis show a near 160-day periodicity in the photospheric 
magnetic flux during solar cycle 21 (but not in solar cycle 

22) and suggest a probable causal relation to the well-known 
Rieger-type quasi-periodicities. Given our data analyses on 
mid-term quasi-periodicities in CMEs and solar X-ray flares, 
it is very likely that during the recent solar maximum (cycle 

23) , the daily total solar magnetic flux of the MWTF, for 
example, should reveal Rieger-type periodicities. The results 
of Ballester et al. (2002) are consistent with the notion that 
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the emergence of subphotospheric magnetic flux may be trig- 
gered or modulated quasi-periodically by large-scale Rossby 
waves trapped in the solar equatorial zone (Lou 2000b). 



6 DISCUSSION AND SUMMARY 

In terms of the interrelation between mid-term quasi- 
periodicities of equatorially trapped Rossby waves and 
Rieger-type periodicities in solar flare related activities, we 
would offer several physical ideas here to address the in- 
teresting questions raised by Ballester et al. (2002) regard- 
ing the mechanism proposed by Lou (2000b). Due to the 
solar differential rotation, subsurface mean magnetic fields 
are wrapped around the solar equator with growing intensi- 
ties. By turbulent convections and MHD instabilities, sub- 
surface magnetic fluxes will emerge in general and sunspots 
or sunspot groups will form in particular. Magnetic activites 
above the photosphere are capable of disturbing the photo- 
sphere and sub-photospheric layers through MHD processes. 

We thus have in mind a "dynamic feedback scenario" 
advanced earlier. Specifically, powerful solar flares frequently 
occurring in magnetic active regions along the two usual 
belts across the solar equator keep stirring the photosphere 
and exciting Rossby waves that in turn either trigger or 
modulate the emergence of subphotospheric magnetic flux. 
To sustain such a "dynamic feedback cycle" , initiations and 
positive feedbacks above a certain energetic "threshold" are 
necessary. This may explain why Rieger-type periodicities 
in the emergence of magnetic flux, sunspot areas, and high- 
energy flares etc. are detected within a few years around the 
solar maxima when solar activity levels are high. This also 
implies possible changes in the most dominant periodicities 
with different physical conditions for different solar maxima. 

Practically, however, it is not easy to estimate this 
"threshold" , because several unknown parameters and pro- 
cesses are involved. Regarding possible values of m and n 
in the solar Rossby wave theory of Lou (2000b), we already 
noted that for ~5 or 6 active regions in the solar disk facing 
us, m would be ~10 or 12 around the circumference. Sim- 
ilarly, equatorially trapped Rossby-type waves with n — 2 
may be more favorably excited given the presence of two 
belts of sunspot groups north and south of the equator. 
In the sense of triggering magnetic avalanches, a sustained 
passage of Rossby-type waves can increase statistically the 
chance of flare occurrence in preexisting vulnerable mag- 
netic complexes and thus impose grossly quasi-periodicities 
to various diagnostics associated with solar flares. 

In summary, by Fourier spectrum and Morlet wavelet 
analyses, we identified quasi-periods in data series of CMEs, 
X-ray solar flares (class>M5.0), and daily averages of Ap in- 
dex during the time span from February 5, 1999 to February 
10, 2003 for the solar maximum of cycle 23. CME periods at 
~ 358±38, ~ 272±26, and ~ 196±13 days appear to corre- 
spond well to Ap index periods at ~ 364±39, ~ 273±26, and 
~ 187±12 days, although the Ap index period of ~ 364±39 
days is somewhat less significant; we interpret these quasi- 
periodic correlations in terms of CME interactions with the 
Earth's magnetosphere. In addition to a significant period 
at ~ 259 ± 24 days, X-ray flare data contains the familiar 
Rieger-type periodicities shorter than ~ 157 ±11 days and 
so forth. Within error ranges, this ~ 259 ± 24-day period in 



X-ray flare data may relate to CME ~ 272 ± 26-day and Ap 
~ 273 ± 26-day periods. We provide estimates for periodic 
timescales of equatorially trapped Rossby-type waves based 
on the work of Lou (2000a, b) and discuss a few qualita- 
tive aspects of a "dynamic feedback cycle" to sustain such 
waves. A coherent picture for solar Rossby waves is yet to 
come given several preliminary results of their circumstan- 
tial detections in recent years (Beck et al. 1998; Kuhn et al. 
2000; Ulrich 2002). More observations and theoretical inves- 
tigations are needed to further understand these mid-term 
quasi-periodicities in solar as well as solar-terrestrial activi- 
ties. 

Finally, we emphasize the strong obvious periodicity of 
187 ± 12-days in daily averaged Ap index that can be visu- 
ally picked out in Fig. 3a. We tentatively identified a CME 
counterpart of period 196 ± 13 days as the cause. A fur- 
ther understanding of this phenomenon is of great interest 
especially in contexts of solar terrestrial interaction. 
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